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We study the properties of spin systems realized by cold polar molecules interacting via dipole- 
dipole interactions in two dimensions. Using a spin wave theory, that allows for the full treatment 
of the characteristic long-distance tail of the dipolar interaction, we find several anomalous features 
in the ground state correlations and the spin wave excitation spectrum, which are absent in their 
counterparts with short-range interaction. The most striking consequence is the existence of true 
long-range order at finite temperature for a two dimensional phase with a broken f/(l) symmetry. 
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The foundation for understanding the behavior and 
properties of quantum matter is based on models with 
short-range interactions. Experimental progress in cool- 
ing polar molecules []| and atomic gases with large mag- 
netic dipole moments [2] has however increased the in- 
terest in systems with strong dipole-dipole interactions. 
While many properties of quantum systems with dipole- 
dipole interactions derive from our understanding of sys- 
tems with short-range interactions, the dipole-dipole in- 
teraction can give rise to phenomena not present in their 
short-range counterparts. Prominent examples are the 
description of dipolar Bose-Einstein condensates, where 
the contribution of the dipolar interaction cannot be in- 
cluded in the s-wave scattering length 3: s-nd the ab- 
sence of a first order phase transition with a jump in 
the density [J|. In this letter, we demonstrate anomalous 
behavior in two dimensional spin systems with dipolar in- 
teractions realized by polar molecules in optical lattices. 

A remarkable property of cold polar molecules con- 
fined into two dimensions is the potential formation of 
a crystalline phase for strong dipole-dipole interactions 
[a, y. In contrast to a Wigner crystal with Coulomb 
interactions 7], the crystalline phase exhibits the con- 
ventional behavior expected for a crystal realized with 
a short-range repulsion and the characteristic Ijr^ be- 
havior of the dipole interaction can be truncated at dis- 
tances involving several interparticle separations. Several 
strongly correlated phases have been predicted, which be- 
have in analogy to systems with interactions extending 
over a finite range, such as a Haldane phase t8|]. super- 
solids [9|, llO| , pair supersolids in bilayer systems ll| , va- 
lence bond solids [l2| , as well as p-wave superfluidity 13 1 , 
and self-assembled structures in multilayer setups 14 1. 



On the other hand, it has recently been demonstrated 
that polar molecules in optical lattices are also suitable 
for emulating quantum phases of two dimensional spin 
models [l5l-[l7|. 

Here, we demonstrate that such spin models with 
dipole-dipole interactions exhibit several anomalous fea- 



tures, which are not present in their short-range coun- 
terparts. The analysis is based on analytical spin wave 
theory, which allows for the full treatment of the l/r^ 
tail of the dipole-dipole interactions. We find that the 
excitation spectrum exhibits anomalous behavior at low 
momenta, which gives rise to unconventional dynamic 
properties of the spin wave excitations. Remarkably, 
we derive from this anomalous behavior the existence of 
a long-range ordered ferromagnetic phase at finite tem- 
peratures; this finding is consistent with the well-known 
Mermin- Wagner theorem as the latter does not exclude 
order for interactions with a 1/r" tail, where a < 4 |18l - 
|20|. Finally, we show that the dipole-dipole interaction 
gives rise to algebraic correlations even in gapped ground 



states, in agreement with recent predictions [21|, [2 




FIG. 1. (a) Mean-field phase diagram for the XXZ model with 
dipolar interactions, where tan 6 is the ratio between the XY 
and the Ising spin couplings, (b) Ground state energy per 
particle: the dashed lines show the mean-field predictions, 
while the solid lines include the contributions from the spin 
waves. At the critical values 6c and 6c, the ground state 
energy exhibits the jump Acc ~ 0.14 J and AEc « 0.06 J, 
indicating the potential formation of an intermediate phase. 

We focus on a set up of polar molecules confined into 
two dimensions in a square lattice, with each lattice site 
filled by one polar molecule. A static electric field applied 
along the z direction splits the rotation levels, and allows 
us to define a spin 1/2 system by selecting two states in 
the rotational manifold. Then, the Hamiltonian reduces 



to a XXZ model with dipole-dipole interaction between 
the spins [16| 

ja3 ^ cos 9 SfS^ + sin elSfS^ + SfSf) 

^=1^5: — ' iR.-R 3^ — -^- w 

Here, the first term accounts for the static dipole-dipole 
interaction between the different rotational levels with 
strength J cos 9 , while the last term describes the virtual 
exchange of a microwave photon between the two polar 
molecules with strength Jsm9, and a denotes the lattice 
spacing. The dependence of the couplings J and 9 on 
the microscopic parameters is discussed in Ref. IQ, ll7, 
12311 and the one-dimensional version of this model has 



recently been studied in Ref. [24 [ . 

Before analyzing this spin model on the square lat- 
tice, we present a summary of the phase diagram for 
its counterpart with nearest neighbor interactions only. 
Then, the phase diagram is highly symmetric and ex- 
hibits four different phases: (i) an Ising antiferromag- 
netic phase (I-AF) for — 7r/4 < 9 < 7r/4 with an excita- 
tion gap, (ii) an XY antiferromagnetic phase (XY-AF) for 
7r/4 < 9 < 37r/4 with a linear excitation spectrum, (iii) 
an Ising ferromagnetic phase (I-F) for 37r/4 < 9 < 5tt/A 
with an excitation gap, and finally (iv) a XY ferromag- 
netic phase (XY-F) for bn/A < 9 < In/A with a linear 
excitation spectrum. 

Next, we analyze the modifications of the phase dia- 
gram due to dipole-dipole interactions between the spins 
within mean-field theory. The main influence is the re- 
duction of the stability for the antiferromagnetic phases, 
as the next nearest neighbor interaction introduces a 
weak frustration to the system. The ground state en- 
ergy per lattice site within mean-field reduces to ei_AF = 
J cos 9 eK/4 and Cxy-af = Jsin^ eK/4 for the antiferro- 
magnetic phases. The summation over the dipole interac- 
tion reduces to a dimensionless parameter ek ~ —2.646, 
which is related to the dipolar dispersion 
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at the corner of the Brillouin zone K = (tt/o, 7r/a). 
In turn, the ferromagnetic phases are enhanced with 
a mean-field energy ei_p = J cos 9 eo/4 and Bxyf = 
Jsin6'eo/4 with eg w 9.033. The modifications to the 
phase diagram are shown in Fig. [T] first, the Heisenberg 
points at 9 — 7r/4, 57r/4 are protected by the SU(2) sym- 
metry and still provide the transition between the Ising 
and the XY phases. However, the transitions from the 
ferromagnetic towards the antiferromagnetic phase are 
shifted to the values 9c — arctan(eK/eo) ~ — O.lvr and 
9c = TT + arctan(eo/eK) ~ 0.67r. 

The dipole dispersion eq in Eq. ([2]) converges very 
slowly due to the characteristic power law decay of the 
dipole-dipole interaction. It is this slow decay, which 



will give rise to several peculiar properties of the system. 
Therefore, we continue first with a detailed discussion of 
this dipolar dispersion. The precise determination of Cq 
is most conveniently performed using an Ewald summa- 
tion [7| , which transforms the summation over the slowly 
converging terms with algebraic decay into a summation 
of exponential factors, i.e.. 



^27ra|q| erfc(a|q|/2\/7r) + An 
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with erfc(x) the complementary error function. The im- 
portant feature of the dipole dispersion is captured by 
the first term in Eq. ([3]), which gives rise to a linear and 
nonanalytic behavior eq ^ eo — 27ra|q| for small values 
q <^ 1/a, while all remaining terms are analytic. It is 
this linear part, which will give rise to several uncon- 
ventional properties of spin systems in 2D with dipolar 
interactions, and is a consequence of the slow decay of 
the dipole-dipole interaction. The summation in the last 
term converges very quickly and guarantees the periodic- 
ity of the dipolar dispersion. The quantitative behavior 
is shown in Fig. [2^, and the numerical efficient determi- 
nation provides eo ~ 9.033, and eK = (l/v^ — l)eo ~ 
-2.646 

Next, we analyze the excitation spectrum above the 
mean-field ground states within a spin wave analysis. 
The spin wave analysis is well established 25|, |26| , and its 
application for a spin system with dipolar interaction is 
straightforward. Details of the calculation are presented 
in the supplementary material for the antiferromagnetic 
XY phase. The results are summarized in Table U and 
shown in Fig. [2l In the following, we present a detailed 
discussion for each of the four ordered phases. 
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FIG. 2. Spin wave excitations with T = (0,0), M = (0,7r/2), 
and K — {n/a, n/a) for different 9 angles, (a) Spectrum of 
the I-F phase which also shows the behavior of the dipolar 
dispersion eq for 9 = — 37r/4, see red line, (b-d) Spectrum for 
the XY-F, I-AF and XY-AF phases. Each red line is a critical 
excitation spectrum indicating an instability. 

Ising ferromagnetic phase: The ferromagnetic mean- 
field ground state is twofold degenerate with all spins 



Ground state a Spin wave excitation spectrum E'^ 



Ground state energy per spin e^ 



I-F 

XY-F 

I-AF 

XY-AF 



3 J cos Qen 1 f da ^ , . J cos Oen 

SJsin^eo 1 f dq ^ 
-^—+2 1^^- 



J (sin Oeq — cos 6eo) 

Jy^sm9{eq — eo){cos9eq — sin^eo) 

J^{sm8eq+K. — cos^eK)(sin^eq — cos6'eK) 

J^sin^(eq+K — eK)(cosS£q — sin^ex) 

TABLE I. Spin wave excitation spectrum E^ and ground state energy ea- 
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— i — + 27 ^^'^ 



either point up or down, and is the exact ground state 
for 9 = n, i.e., jG) = fl^ ||)j. Within the spin wave anal- 
ysis, the ground state is not modified and the excitation 
spectrum reduces to E'^^ , see Table HI The spin waves 
exhibit an excitation gap A: (i) approaching the Heisen- 
berg point at 6* = — 37r/4, the excitation gap vanishes, 
indicating the instability towards the XY ferromagnet, 
(ii) in turn, for antiferromagnetic XY couplings, the gap 
is minimal at K, vanishes at the mean-field transition 
point 9c and drives an instability towards the formation 
of antiferromagnetic ordering. 

In contrast to any short-range ferromagnetic spin 
model, the dispersion relation E^^ is not quadratic for 
small momenta, but rather exhibits a linear behavior, 
i.e., E'^^ ^ i^Q^ + fi.c|q| with velocity c — —2TraJsm9/h, 
which is a consequence of the dipolar interaction in the 
system. This anomalous behavior strongly influences the 
dynamics of the spin waves. The dynamical behavior of 
a single localized spin excitation is shown in Fig. |3^ for a 
Gaussian initial state. In order to probe the linear part in 
the dispersion relation, the width a of the localization is 
much larger than the lattice spacing a, and therefore, the 
dynamics is well described by a continuum description. 
Instead of the conventional quantum mechanical spread- 
ing, one finds a ballistic expansion of a cylindrical wave 
packet with velocity c. In addition, the dipole-dipole 
interaction also strongly influences the correlation func- 
tion. Within conventional perturbation theory, we find 



algebraic correlations {SfS^ 



1/ 



This algebraic 



decay of correlations even in gapped systems is a pecu- 
liar property of spin models with long-range interactions 



21,^ 



XY-ferromagnetic phase: Here, the spins are aligned in 
the xy plane. Within the spin wave analysis, we obtain 
the excitation spectrum E^'^'^ and the modified ground 
state energy Cxyf- In the low momentum regime, the 
dispersion relation behaves as E^^'' ~ \/|qf, in contrast 
to the well known linear Goldstone modes for the broken 
C/(l) symmetry. This anomalous behavior is a peculiar 
property of the dipolar interaction, and the most cru- 
cial consequence is the existence of long-range order for 
the continuous broken symmetry at finite temperatures 
even in two dimensions (l9|. This property follows im- 
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FIG. 3. Time evolution for localized spin excitations de- 
scribed by the Gaussian wave packet i/)o(r) = e"'''' ""^ /Vtvo^ 
with cr 3> a in the continuum description, (a) For a linear 
dispersion cjqj in the I-F phase, the dynamics is described 
by cylindrical symmetric wave packets (see inset) traveling 
with velocity c, instead of the conventional quantum mechan- 
ical spreading for massive systems, (b) For an anomalous 
dispersion with Qy^qJ in the XY-F phase, the behavior at 
long times t ^ y/^Q: reduces to a scaling function ^{z) via 
|V'(a;,''")|'^ ~ ^{x/t ^ 1/2)/''"^ (see inset) using rescaled time 
T = ta/y/a and space x = \r\/a coordinates. It describes a 
cylindrical symmetric wave front with velocity a^/a. 



mediately from the above spin wave analysis: the order 
parameter reduces to to = Am — 1/2 = {Sf)/h, where 
Am accounts for the suppression of the order parameter 
by quantum fluctuations. Within spin wave theory, it 
reduces to (Ato, — {alai)) 



Am— I — 



cos oer. 



smi 



4.Er. 



^^^^cothr^Vi 



2r 



This expression is finite and small: at T := 0, the inte- 
grand behaves as ~ l/-^/|q| and we find a suppression of 
the order Am w 0.008 at 6* = —tt/2. The smallness of this 
corrections due to quantum fluctuations is a good justifi- 
cation for the validity of the spin wave analysis. On the 
other hand, even at finite temperatures, the low momen- 
tum behavior of the integrand takes the form ~ 2^/|q|, 
and provides a finite contribution in contrast to a con- 
ventional Goldstone mode, which provides a logarithmic 
divergence. 

The appearance of a long-range order at a finite tem- 
perature for a ground state with a broken U{1) sym- 
metry is a peculiar feature of dipole-dipole interactions, 
which renders the system more mean-field like. The sys- 



correlation function 



r = o 



< T < Te 



Tc<T 



(Sf Sy + SfSJ) - m^ 



-5/2 
-3/2 



TABLE II. Correlation functions in the XY-F phase predicted 
by the spin wave analysis and high temperature expansion. 



tern therefore exhibits a finite temperature transition at 
a critical temperature Tc into a disordered phase; such a 
behavior is consistent with the classical XY model with 
dipolar interactions [19[ . The correlation functions deter- 
mined within spin wave theory and a high temperature 
expansion are summarized in Table HIl Note, that the spin 
wave analysis neglects the influence of vortices. This is 
well justified here, as the dipolar interactions give rise 
to a confining of vortices, i.e., the interaction potential 
between a vortex-antivortex pair increases linearly with 
the separation between the vortices. 

The spin wave dynamics caused by the anomalous dis- 
persion relation ^ vlq| is shown in Fig.[3jD for a Gaussian 
wave packet of width a. Interestingly, the propagation 
velocity of the wave packets is proportional to y^ and 
thus faster for broad wave packets, in contrast to the 
usual dispersion dynamics. This is a consequence of the 
group velocity Vq ^ l/y^q| which is large for the small 
momentum components involved in the broad wave pack- 
ets. 

Ising antiferromagnetic phase: Next, we focus on the 
antiferromagnetic phases and start with the I-AF ground 
state. Again, the ground state is twofold degenerate 
on bipartite lattices. We choose the ground state with 
spin up on sublattice A and spin down on sublattice B, 
i.e., \G) = UieA\'t)^I\j<^B\i)j■ The spin wave analy- 
sis is straightforward (see supplementary information), 
and we obtain the spin wave excitation spectrum E*'' 
and ground state energy ei_AF, see Table HI The sys- 
tem exhibits an excitation gap as expected for a system 
with a broken Z2 symmetry. However, the dipole in- 
teractions give rise to an anomalous behavior at small 
momenta similar to the ferromagnetic Ising phase with 
^i^-AF_^j-AF ~ — sin6'|q|. Consequently, the dynamics of 
spin waves at low momenta is analogous to the Ising fer- 
romagnet, see Fig.[3l Within spin wave theory, we obtain 
that the antiferromagnetic correlations ((— l)'^-'S'f 5^ 
and the ferromagnetic correlations {S^ S^) decay with 
the power law ^ l/|rp with /3 — x,y,z determined 
by the characteristic behavior of the dipole-dipole in- 
teraction. The excitation gap vanishes approaching the 
mean-field critical point 9c towards the XY- ferromag- 
netic phase, and also while approaching the Heisenberg 
point at 6* = tt/A. For the latter, the qualitative behavior 
of the excitation spectrum changes drastically within a 
very narrow range of 0, see Fig. [2j;. 

XY antiferromagnetic phase: Finally, we analyze the 
properties of the antiferromagnetic XY phase. In con- 



trast to the ferromagnetic XY phase, the excitation spec- 
trum ET^'^^ exhibits the conventional linear Goldstone 
mode, see Fig. [2ji. This can be understood, as the an- 
tiferromagnetic ordering introduces a cancellation of the 
dipolar interactions, and provides a behavior in analogy 
to its short-range counter part: true long-range order 
exists only at T = 0, while at finite temperature the 
system exhibits quasi long-range order and eventually 
undergoes a Kosterlitz-Thouless transition for increasing 
temperature. Nevertheless, the dipole-dipole interactions 
give rise to the characteristic algebraic correlations, e.g., 
{{—ly^^ S^Sp ^ l/|r|^ for the antiferromagnetic trans- 
verse spin correlation at zero temperature. 

Finally, we comment on the transitions between the 
different phases. The spin wave analysis predicts, that 
the excitation spectrum for each phase becomes unsta- 
ble at the mean-field critical points: For the Heisenberg 
points a.t 6 — 7r/4, 57r/4, such a behavior is expected due 
to the enhanced symmetry and one indeed finds, that at 
the critical point, the excitation spectrum from the Ising 
phase coincides with the spectrum from the XY ground 
state. Consequently, the spin waves provide the same 
contribution to the ground state energy, see Fig. [T]d. In 
turn, at the instability points Oc and 9c-, the excitation 
spectrum from the antiferromagnetic phase is different 
from the spectrum for the ferromagnetic F phase. Con- 
sequently, the ground state energy within the spin wave 
analysis exhibits a jump, see Fig. [1^. Such a behav- 
ior is an indication for the appearance of an intermediate 
phase. However, this question cannot be conclusively an- 
swered within the presented analysis due to the limited 
validity of spin wave theory close to the transition points. 
However, the appearance of a first order phase transition 
can be excluded by arguments similar to Ref. [J|. 
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SUPPLEMENTARY INFORMATION 

We present the derivation of the spin wave excitation 
spectrum within the spin wave analysis. The basic ap- 
proach is to start with the ground states exhibiting per- 
fect order, which are the correct ground states for the 
classical model at the four points 6 — 0, ±7r/2,7r. Then, 
we introduce bosonic creation and annihilation operators 
creating a spin excitation above the ground state accord- 
ing to the Holstein-PrimakofT transformation. The spin 
Haniiltonian then reduces to a Bose-Hubbard model. In 
lowest order, we can ignore the interactions between the 
bosonic particles, and obtain a quadratic Hamiltonian 
in the bosonic operators, which is diagonalized using a 
Bogoliubov-Valantin transformation. The latter trans- 
formation deforms the ground state and introduces fluc- 
tuations into the system. 

In the following, we demonstrate the spin wave analy- 
sis for the most revealing case: the antiferromagnetic XY 
phase. The generalization to the other ground states is 
straightforward. Without loss of generality, we choose 
the antiferromagnetic order to point along the x di- 
rection. The square lattice is bipartite, and we de- 
note the two sublattices as A and B. Then, the anti- 
ferromagnetic mean- field ground-state is given by \G) = 
Y\i^A \^iY\ii^B \~^)j with spins on sublattice A point- 
ing in the negative x direction, i.e., S^ |<— ) = — fi,/2 |<— ), 
and spins on sublattice B pointing in the positive x di- 
rection (5"^ j-^-) = h/2 |— >■)). Excitations on sublattice A 
are created by flipping a spin with the ladder operator 
S^^ — S^ — iS"^, while excitations on sublattice B are cre- 
ated via S^~ = S^ -\-iSy . We apply a Holstein-Primakoff 
transformation to bosonic operators 



2K 



](f{ni 



^r-^K 



a\)e^^ 



' ^{nt 



where the phase e*^^' 



*^^' accounts for the 



sublattice-dependent sign with K = (7r/a,7r/a). The fac- 
tor ^{rii) = 1 — ni is introduced to guarantee bosonic 
commutation relations for the operators a^. Here, we are 
interested in the leading order of the spin wave expansion, 
and can therefore set ip{ni) « 1. The bosonic operators 
reduce to 

and the number operator rii — a\a^ — 2 ^ SxC^^^' /H. 

Expanding the spin Hamiltonian in terms of the 
bosonic operators leads to a Bose-Hubbard Hamiltonian 
for the spin wave excitations. In leading order, we can 
neglect the interactions between the bosons and obtain 



H 
~J 



4e 



6N 1 V- r t t 



(4) 



with R, 



Ri — Rj, N the number of lattice sites, 



and the couphng the terms Xij = cos6' + sin6'e*^^'J 
and Tjij = cos 9 — sin Oe'^^^'^' including the antiferromag- 
netic ordering. Introducing the Fourier representation 
tti — X]q^q6~*'*^'/v^, the terms involving the bosonic 
operators in Eq. Q reduce to 

- ^ (cos 9€q + sin 6 tq+K - 2 sin 06^) [al^a^ + a^a^^ 



+ (cos Eq - sin 6 Eq+K) ( ftqft-q + Aqfl-, 



The diagonalization of this Hamiltonian is straightfor- 
ward using a standard Bogoliubov transformation with 



form 



q^-q 



Then, the Hamiltonian takes the 



H 



3 JN sine EK 



\ TT-XY-AF 



6^&q 



(5) 



with the spin- wave excitation spectrum E^^'^^. In addi- 
tion, the coefficients for the Bogoliubov transformation 
are given by 




cos 9eq + sin 9 (eq+K - Sek) 



±1 (6) 



with £„ 



rpXY-AF 



I J . The property «,; 



1 asserts 



that the transformation is canonical. In addition, the 
ground state obeys the condition 6q |vac) = 0, and the 
ground state energy per spin at zero temperature T = 
reduces to Cxy-af, see Table HI 

We are now able to check the validity of the spin 
wave approach self-consistently: the deformation of the 
ground state by the spin wave analysis provides a sup- 
pression of the antiferromagnetic order m = Attt, — ij = 
(S'fe*^^')/S, and thus 



'I'D J vq 



+ {2vl + l)/q , (7) 



where /q = ih'^^^ = [exp(£;^^-*^/r) - l] accounts 
for the thermal occupation of the spin waves. At zero 



temperature T = 0, this expression converges and we 
obtain Am sa 0.03 for 6* = 6'c as well as Am = 0.39 for 
sa ^. In turn, at finite temperatures T > 0, the low 
momentum behavior of the integrand scales as |q|~^, and 
therefore Am. diverges logarithmically: the long range 
order is destroyed by the thermal spin wave fluctuations, 
and gives rise to the well-known quasi long-range order 
in analogy to short-range XY models. 

Finally, the spin wave analysis also allows us to ana- 
lyze the correlation functions Caa^^ij) = (Sf S'"e*^^'j). 
Using the translational invariance of our system, the cor- 
relation functions reduce to 



.(r) 






Caaiq + 'K) 



-zqr 



(8) 



with CQQ(q) = {SaS'^ ). Combining the Holstein Pri- 
makoff transformation Eq. @ and the Bogoliubov trans- 
formation the correlation functions can be expanded in 
terms of the coefficients Uc, and «„ 



1 f J£ 

C^^(q+K) = -(Wq+K+l^q+K)^COth ( -^ 
Cyy{c\ + K) = -(Wq - V^f COth ( -^ 



(9) 



The long distance behavior |r| — )■ oo of the correlation 
function is determined by the low momentum behavior 
of the above expressions 



(Wq+K 



f'q+Kj 



(Wq 



|q| + const. 
1 



(10) 



and describes the leading non-analytic part. The latter 
can be replaced using the following relation, which de- 
rives via an Ewald summation. 



oiqRj 






(11) 



for |q| -!■ and 7 > -2; (for 7 = 0,2,4,... the left side 
is replaced by jqj''' log |q|). At zero temperature T = 0, 
the integration in Eq. (|8]) is straightforward and provides 
the scaling behavior Czz{'r) ^ [rj^'^ and Cj^j^(r) ~ |r|~^- 



